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Abstract. The Euler-Poisson system is a fundamental two-fluid model to de- 
scribe the dynamics of the plasma consisting of compressible electrons and a 
uniform ion background. By using the dispersive Klein-Gordon effect, Guo 
[5] first constructed a global smooth irrotational solution in the three dimen- 
sional case. It has been conjectured that same results should hold in the two- 
dimensional case. The main difficulty in 2D comes from the slow dispersion 
of the linear flow and certain nonlocal resonant obstructions in the nonlincar- 
ity. In this paper we develop a new method to overcome these difficulties and 
construct smooth global solutions for the 2D Euler-Poisson system. 

1. Introduction 

The Euler-Poisson system is a fundamental two-fluid model used to describe the 
dynamics of a plasma consisting of moving electrons and ions. Since the heavy 
ions move much more slowly than the light electrons, one can treat them as two 
independent compressible fluids which only interact through their self-consistent 
electromagnetic fields. In the simplest approximation, the ions are immobile and 
uniformly distributed in space, providing only a background of positive charge for 
the electrons. Neglecting magnetic effects, the equations of motion describing the 
dynamics of the compressible electron fluid is then given by the following Euler- 
Poisson system in {t,x) G [0,oo) x M.'^, 

' dtn + V ■ (nu) = 0, 

u + (u • V)u) -I- Vp( i) = e7xV0, (1.1) 
A(/) = 4:Tre{n — no). 

Here n = n{t, x) and u = u(t, x) denote the density and average velocities of 
the electrons respectively. The dimension d will mostly be either d — 3 ot d = 2 
and we shall state clearly the dimensional dependence of the results. The state 
of the plasma is completely specified once the density and average velocities are 
known. It is useful to keep in mind that such a fluid model description is only 
accurate when the plasma velocity distribution is close to the Maxwell-Boltzmann 
distribution so that the microscopic details of each individual particles are averaged 
out. The symbol e and nie denote the unit charge and mass of electrons. The first 
equation in (jl.ip is simply the mass conservation law, while the second equation 
comes from the momentum balance. The pressure term p{n) is assumed to obey 
the polytropic 7-law, i.e. 

p{n) = An'', (1.2) 

where A is a constant and 7 > 1 is usually called the adiabatic index. The last 
equation in (|1.1|) is the Poisson equation (or Gauss law) which computes the electric 
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field self-consistently through the charge distribution. We assume here that at the 
equihbrium the density of ions and electrons are both a constant denoted by ng. 
The dependence on ions only enter through the term 47re(ri, — uq) which is the net 
charge generating the electric field. It is in this sense that the system (jl.ip is the 
simplest two-fluid model for plasmas since the ions arc merely treated as constant 
uniform sources. The natural boundary condition for the electric potential is a 
decaying condition at infinity, i.e. 

lim (f>{t,x) = 0. (1.3) 

|a;|— ^oo 

Throughout the rest of this paper, we shall consider an irrotational flo"v\3, i.e. 

Vxu = 0. (1.4) 

By a simple computation, it is easy to check that the irrotational condition is pre- 
served for all time. The system (|l.ip - (|1.4p describes the dynamics of an irrotational 
compressible Eulerian electronic fluid moving against the ionic forces. 

Our main objective is to construct global in time smooth solutions for the Euler- 
Poisson system. This is not an easy task since this is a hyberbolic conservation law 
with zero dissipation which has remain largely unsolved and no general theory is 
available. The Euler-Poisson system has close connection with the well-known com- 
pressible Euler equations. Indeed in (jl.ip if the electric field term Vcj) is dropped, 
one recovers the usual Euler equations for compressible fluids. In |14j . Sideris con- 
sidered the 3D compressible Euler equation for a classical polytropic ideal gas with 
adiabatic index 7 > 1. For a class of initial data which coincide with a constant 
state outside a ball, he proved that the lifespan of the corresponding solu- 
tion must be finite. In [12] Rammaha extended this result to the 2D case. For 
the Euler-Poisson system, Guo and Tahvildar-Zadeh [6] established a "Siderian" 
blowup result for spherically symmetric initial data. Recently Chae and Tadmor 
[2] proved finite-time blow-up for solutions of a class of pressureless attractive 
Euler-Poisson equations in M", n > 1. These negative results showed the abundance 
of shock waves for large solutions. 

The point of departure in the analysis of (|1.1|) from that of the compressible 
Euler equations begins with understanding the small irrotational perturbations of 
the equilibrium state n = hq, u = 0. For the 3D compressible Euler equation with 
irrotational initial data (n, u) {ep + no, ev), where p e 5(M^), v e 5(M^)^ (5(M^) 
is the usual Schwartz space), Sideris [15] proved that the lifespan of the classical 
solution > exp(C/e). For the upper bound it follows from his previous paper 
[14j that Te < exp(C/e^) under some mild conditions on the initial data. Sharper 
results are also available. For initial data which is spherically symmetric and is 
smooth compact e-perturbation of the constant state, Godin 0] obtained by using 
a suitable approximation solution the precise asymptotic of the lifespan as 

lime log = T*, 

where T* is a constant. These results rely crucially on the observation that after 
some simple manipulations, the compressible Euler equation in rescaled variables 
is given by a vectorial nonlinear wave equation with pure quadratic nonlinearities. 
The linear part of the wave equation decays at most at the speed which in 

^In the two dimensional case, one can regard in the usual fashion u = {ui,U2,0), for which 
the condition V X u = simplifies to dx^U2 — dx^'i^i- 
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3D case is not integrable. Unless some additional structure e.g. the null condition 
[U [9] is imposed on the nonlinearity, one cannot in general expect global existence 
of small solutions. For the Eulcr-Poisson system the situation for small irrota- 

tional perturbations of the equilibrium state is quite different. By a straightforward 
computation, the linear part of the Euler-Poisson system for irrotational flows has 
the form 

^ttn - ;^p'(ri,o)An + uj'^n = 0, 
dttu - :^p'{no)Au + ujIvl = 0. 



Here ujp = ^ Aire^nQ/nig^ is so called plasma frequency. This new term uj'^ which is 
absent in the pure Euler case makes the linear part of the Eulcr-Poisson system into 
a Klein-Gordon system for which the linear solutions have an enhanced decay of 
{l+t)-'^/'^. This is in sharp contrast with the pure Euler case for which the decay 
is only ^-(''-i)/^. Note that in d = 3, (1 t)-'^/'^ = (1 + t)"^/^ which is integrable 
in t. In a remarkable paper [5], by exploiting the crucial decay property of the 
Klein-Gordon flow in 3D, Guo [S] modified Shatah's normal form method [13] and 
constructed a smooth irrotational global solution to (jl.ip around the equilibrium 
state (riQ, 0) for which the perturbations decay at a rate Cp- (l+t)"^ for any 1 < p < 
3/2 (here Cp denotes a constant depending on the parameter p). Note in particular 
that the sharp decay t~^^'^ is marginally missed here due to a technical complication, 
namely in employing one of the — )■ L°° decay estimates, the nonlinear term has a 
Riesz-type singular operator in the front and one is forced to take p > 1 in bounding 
the corresponding part. 

As was already mentioned, the purpose of this paper is to construct smooth 
global solutions to for the two-dimensional case. One can already sense a bit 
that this situation is more challenging than the d ~ Z case since the linear solutions 
to the Klein-Gordon system in d = 2 decays only at (1 -I- t)~^ which is not inte- 
grable, in particular making the strategy in [5] difficult to apply. As we shall see 
shortly, the Euler-Poisson system after writing in terms of the perturbed variables 
can be recast into a quasi-linear Klein-Gordon system with nonlocal nonlineari- 
ties. To understand and spell out the main difficulties associated with the analysis 
of this nonlocal quasi-linear Klein-Gordon system, wc first review some known 
results on nonlinear Klein-Gordon equations. For general scalar quasi-linear Klein- 
Gordon equations in 3D with quadratic type nonlincarities, global small smooth 
solutions were constructed independently by Klainerman [8] using the invariant 
vector field method and Shatah [13j using a normal form method. There are es- 
sential technical difficulties to employ Klainerman's invariant vector field method 
to the Klein-Gordon system associated with (|l.ip . Namely after reformulating in 
terms of perturbed the variables, the nonlinear part of the system (jl.ip has a Riesz 
type nonlocal term which comes from solving the Poisson equation for the electric 
field in i.e. a term of the form 

{—/S)^'^dxidxj (quadratic terms). 

The Klainerman invariant vector fields consist of infinitesimal generators which 
commute well with the linear operator 9tt — A + 1. The most problematic part 
comes from the Lorentz boost fioj = tdx, + Xjdt- Whilst the first part tdx^ com- 
mutes naturally with the Riesz operator i?y = [~A)~^dxidxj , the second part xjdt 
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interacts rather badly with Rij and produces a commutator which scales as 

[xJ^uR^,]^^t\V\-\ 

After repeated commutation of these operators one obtain in general terms of the 
form |V|~^ which makes the low frequency part of the solution out of control. It 
is for this reason that in 3D case Guo [S| adopted Shatah's method of normal form 
in L'P {p > 1) setting for which the Riesz term i?y causes no trouble. We turn now 
to the 2D Klein-Gordon equations with pure quadratic nonlinearities. In this case, 
direct applications of either Klainerman's invariant vector field method or Shatah's 
normal form method are not possible since the linear solutions only decay at a 
speed of {l+t)~^ which is not integrable and makes the quadratic nonlinearity quite 
resonant. In |16| . Simon and Taflin constructed wave operators for the 2D semilinear 
Klein-Gordon system with quadratic nonlinearities. In [11) . Ozawa, Tsutaya and 
Tsutsumi considered the Cauchy problem and constructed smooth global solutions 
by first transforming the quadratic nonlinearity into a cubic one using Shatah's 
normal form method and then applying Klainerman's invariant vector field method 
to obtain decay of intermediate norms. Due to the nonlocal complication with the 
Lorentz boost which we explained earlier, this approach seems difficult to apply in 
the 2D Euler-Poisson system. 

Our approach in this paper is inspired by the recent work of Gustafson, Nakanishi 
and Tsai [3] on the Gross-Pitaevskii equation of the form 

idt^p = -Ai;+i\ij\^ 

where ip : ^ C is solved with the boundary condition 

\ip{t,x)\ — > 1, as \x\ — )- oo. 

The main objective in is to investigate large time behavior of solutions ip ~ 
small". It is then natural to look at the perturbation u = ip — 1 which satisfies 
the equation 

idtu + Au-2Re{u) ^ F{u), F{u) := + 2\u\'^ + \u\'^u. (1.5) 

To make the problem complex linear, introduce change of variable 

w = /Cu := v/(-A)-i(2 - A)Re{u) + ilm{u). 

Then for v the equation takes the form 

idtv - ^J{~A){2-A)v = -iK.iF{JC-^v). (1.6) 

In [3], Gustafson, Nakanishi and Tsai considered the final data problem for (|1.6p 
(i.e. Cauchy data at t = oo), 

/oo 

where 4> is the final data. Under some decay and smallncss assumptions on (p, they 
constructed a unique global solution for (|1.6p in dimensions d = 2, 3. In particular 
for d = 2, the solutions has a decay 

\\v{t) - e-V-A(2-A)t^||^^ < (1.7) 
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for any e G (0,1) and Ce is a constant depending on e. Note that this marginally 
misses the sharp decay t^^ due to some technical issues which we explain now. The 
first step in [2j is to estimate the first nonlinear iterate, i.e. 

where 

i;i(t)=e-V-A(2-A)t^ 

is simply the linear solution. By a careful space-time phase estimate (see Lemma 
4.1 in [3]), the authors showed that 

here we suppressed the dependence on the final data (f>. The logarithm part of 
t comes mainly from certainly degeneracies of the phase near ^ = 0. Since the 
nonlinearity is quadratic, one immediately observe that due to this logarithm loss 
we cannot close the simple estimates for further iterates (in fact there also some 
problems with low frequencies which we shall not dwell further here). To remedy 
this problem, Gusatfson, Nakanishi and Tsai then introduced a novel normal form 
for the original equation written in the function u (see p.Sp ) 

w = u + P<i{\u\y2), 

where P<i is the usual Littlcwood-Paley projector localized to frequency |^| < 1. 
After this transformation and using Strichartz estimates together with a fixed point 
argument, the authors were able to settle for this logarithmic loss and obtain the 
existence of global solutions, in particular establishing (jl.7p . 

Our analysis of the Euler-Poisson system starts with sharpening the analysis of 
the Gross-Pitaevskii equation in 2D. In [TU], the second author of this paper ob- 
served that provided with some suitable transform and choosing a good working 
space, one can remove the epsilon in (jl.7|) . in particular obtaining the sharp decay 
1/t in 2D. A key observation in this work is that one should obtain the critical 1/t 
decay in the energy estimates of the first nonlinear estimate, and after that only 
energy estimates are needed (in particular, no Strichartz estimates are needed!). 
This intuition is further arrested in our subsequent work [7], where we revisited 
the classical semilinear Klein- Gordon equations in 2D with quadratic nonlinearity 
and also semilinear wave equations in 3D with a quadratic nonlinearity satisfying 
the null condition. Note that in both cases the linear solution decay at the crit- 
ical speed 1/t for which the quadratic nonlinearity is resonant. In that work we 
constructed wave operators for the final data problem and obtain the sharp decay 
of the nonlinear part of the solution. As a byproduct of our analysis in [7], one 
can reinforce the notion of the null condition as an annihilation condition on the 
quadratic interaction of linear solutions. 

In this work, we take the point of view that the Euler-Poisson system is a quasi- 
linear Klein- Gordon system with nonlocal quadratic nonlinearities. We shall con- 
struct smooth global solutions to this system by solving the final data problem 
(i.e constructing the wave operators). As is well-known with quasilinear systems, 
the Strichartz estimates will suffer loss of derivatives unless one considers energy 
estimates. We shall resolve this difficulty by performing a delicate localization in 
time argument. This localization is only carried out for the part of the solution 
which has no resonances and only for the corresponding quasi-linear part of the 
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nonlinearity. We then introduce a new iteration scheme to construct the solution 
as a uniform limit of iterates. The crucial point in the iteration step is to obtain a 
priori uniform sharp decay estimates of i/'" norms of the nonrcsonant part of the 
solution. This is done by a delicate bootstrap argument. Wc should stress that our 
new iteration scheme is quite robust and should have applications to many other 
systems for which the nonlinearity is resonant. To better illustrate our technique, 
we summarize below the main steps of the proof. 

Step 1. Reformulation. For simplicity set all physical constants e, TOc, 47r and 
A to be one. To simplify the presentation, we also set 7 = 3. Define the rescaled 
functions 

n{t/co,x)-no 

u[t,x) = , 

no 

v{t,x) = —u{t/co,x), 
Co 

where the sound speed is co = a/S't-o- For simplicity we shall set no = 1/3 so that 
the characteristic wave speed is unity. Then in terms of the new functions (u, v), 
the Eulcr-Poisson system (|l.ip - (|1.4p takes the form 



(□ + l)u = V • (^(v • V)v + uVuj - dt[v -Vu + uV ■ vj 

(□ + l)v = V('(v • V)m + itV • v) - 9tf (v • V)v + uVu) - VA'^V • (uv), 



or in a more compact notation 

+ = Gi(u,v) 
[(□ + l)v = G2(u,v). 

Step 2. Remove resonances and fine decomposition of the solution. Define 
linear solutions according to p.ip - p.2p . We then identify the resonant part of the 
solution corresponding to the quadratic interaction of the linear solutions. This is 
done by solving the system (see p.5p - p.6p ) 

f(n + l)$i =Gi(/ii,/i2), 

j(n + l)$2 =G2(/ll,/l2) 

with a decaying condition at t = oo. We then decompose the solution as 

ju = u + hi + $1, 

[v = V + /l2 + $2. 

The new variable {u, v) no longer contains resonances and satisfy 

= Gi(?2,v,/i,$), 




= G2(w,v,/i,$), 

where Gi, G2 are corresponding new nonlinearitics. 

Step 3. Sharp decay estimates for hi, /12 and the resonances $1, $2- In this 
step we perform a careful space-time analysis of the bilinear integral for $1, $2- 
Note that the L°° norms of the source terms hi, /12 decay as 1/t which is a standard 

— L°° estimate. The crucial point in our analysis is to establish sharp 1/t decay 
of i/™ norms of $1, $2- This is needed later for the iteration step. 
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Step 4. Localization in time and definition of the iteration scheme. Let x S 
C^iR) be such that x{t) = 1 for < t < 1 and x(i) = for t > 2. For integers 
n > 0, define 

Xnit) = x(2-"t). 

Observe that Xn is locahzed to the regime t < 2". Let u'^^^ = 0, v^"^ = 0. For n > 1 
inductively define 

({a + = x„W(Gn(il("), •) + Gi2(v("), •)) + Gi3(w("-i\ v("-i), /i, $), 

[(□ + l)v(«) = XnW(G2l(ll("\ •) + G22(v("), •)) + G23(w("-1\ v("'l), /l, $). 

Here Gn, G12 denote the quasi-linear part of the nonlinearity in Gi. Similar 
convention applies to G2 (see p.lOp - p.lip for more details). Note that the time 
localization is only applied in the quasi-linear part of the nonlinearity containing 
only the interactions of the non-resonant part of the solution. 

Step 5. Solvability of the iteration system and non-uniform energy estimates. 
The solvability of the iteration system is not immediately obvious, due to the quasi- 
linear nature of the problem. In this respect our strategy is to take advantage of 
the smooth time cutoff function Xn{t) which vanishes for t > 2"+^. In the large 
time regime t > 2"+^, the nonlinearity of the iteration system does not contain 
(?2("),v(")) and we can solve for all t > 2"^^^ by using the fundamental solution 
with source term depending only on the previous iteration (m^""^-*, v^""^-'). After 
this is done we solve backwards the same system for t < 2"+^ with initial data 
(-u("^(2"),v'^"'(2")). We then perform an energy estimate. Denote 

+ ll(v''^W)'lk", Vfc>o. 

Our inductive hypothesis is that 

A„-i(i)<-^, yt>o. 



Here (t) = vT+t^ is the usual Japanese bracket notation. We then use a crude 
energy estimate to obtain 

. , , Kne 

An{t) < 



where Kn is a constant depending on the iteration n. This dependence on the iter- 
ation number n comes from the part where we solve the iteration system backwards 
for t < 2"+i. 

Step 6. Bootstrap and contraction. In this step we upgrade the estimate on 
An and remove the dependence on n. This is done by an infinite-time Gronwall 
argument (see Lemma 16. 1|) which rests on the intuition that nonlinear decay is 
dictated by the sharp decay of the source terms. After obtaining the uniform energy 
estimate and closing the induction hypothesis, we perform the usual contraction 
argument. There is one more small twist in the argument, namely we have to first 
show the strong contraction with a weak time decay but this is enough to 

extract the limiting solution since the true decay {t)~^ from the iterates can be 
passed to the limit for any finite t. 
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Organization of the paper. In section [2] we give the reformulation of the 2D 
Euler-Poisson system and state the main results. In section [3] we extract the non- 
resonant part of the solution and introduce the iteration scheme. Section 2] is 
devoted to the estimates of the resonant part of the solution. In section [5] we prove 
solvability of the iteration system and establish some non-uniform energy decay 
estimates. In section [6] we upgrade the decay estimate of the iterates and establish 
uniform in time energy estimates. In the last section [7l we perform a contraction 
argument and prove the existence and uniqueness of the desired classical solution. 
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2. Reformulation and main results 

We first transform the system (jl.ip - ()1.4p in terms of certain perturbed variables. 

For simplicity set all physical constants e, me, 47r and A to be one. To simplify the 

presentation, we also set 7 = 3 although other cases of 7 can be easily included in 

the discussion. Define the rescaled functions 

n{t/cQ,x)~nQ 

u{t,x) = , 

no 

v{t,x) = —u{t/co,x), 
Co 

■(/'(t, x) ~ (f){t/co, x), 

where the sound speed is cq — VSuq- For simplicity we shall set no = 1/3 so that 
the characteristic wave speed is unity. The Euler-poisson system - (|1.4p in new 
variables take the form 

' (9f u + V • V + V • Vu + uV • V = 0, 
dtv + Vu + (v • V)v + uWu = 3VV', 
dtS/tlj = -iv - iVA^iy • (itv), 

Since the fiow is irrotational V x v = 0, we can then eliminate the electric field 
\/ip using the Poisson equation. Take one more derivative and after some simple 
computation, we arrive at the klein-Gordon form of the 2D Euler-Poisson system 
written for the perturbed variables (u,v): 



(□ + l)u = V • (^(v • V)v + uVuj - dt[v -Vu + uV ■ 

{D + l)v = V(^(v • V)u-(-mV • - dt(^(v ■ V)v + uVuj - VA^^V • (uv). 

(2.1) 

To simplify the discussion later, it is useful to denote the system as 

('(□ + l)« = GiKv) 

|(n + i)v = G2(u,v). ^ ■ ' 
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In Duhamel formulation, we write 



*sin(t-s)(V) /Gi(ii,v)(s) 



ds. 



Note that 



cos(t - T)(V) = cosi(V) cosT(V) + sint(V) sinT(V) 
sm{t-T){V) =sint(V) cosr(V) - cost(V) sinT(V). 



Therefore 



e„.<V,(„..(V, -^ilip (:<?))) ,,3, 

r(v),i„r(v) f!<?;n +»sT(v) f!'<?:m (2,4) 



Consider the linear flow generated by the homogeneous equation 

({a + i)u = o 
|(n + i)v = o 



with initial data 



v(0)j Ui' U(0)j I VI 



It is not difBcult to check that 



u{t)\_( cosi(V) ^ 



ut{t)) l-(V)sint(V) cost{W) I \ui 



Uq 



and similarly 



/ vo 



vVt(i)y l-(V)sinf(V) cost(V) / V^i 
Write the linear flow propagator 

y-(V)sint(V) cost(V)y 

For given data (/i,gi), (/2,52), the wave operator problem for (|2.ip amounts to 
finding the nonlinear solution {u{t), v(t)) such that 

^(-r) (^^^J^) ^ f^'l ' as T ^ oo. (2.7) 
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The convergence here takes place in some Sobolev space if™ which will become 
clear later. Plugging ^3^, dfj]) into ([231), we get 



The system p.Sp will be our main object of study. Here is our theorem. 

Theorem 2.1. Let m > be an integer. Let (fi,f2), (gi, 52 J given functions 
such that 

sup (0'"+'^ E 1^? + E 1^? < oo, (2.9) 

S^"*^ |q|<2 |a|<2 

i=l,2 i=l,2 

where f^, g^ are their Fourier transforms. Let 



Then there exists eg > 0, such that for any < e < eq, there exists a unique global 
solution to (|2.8p written in the form 



^At)J ' ' \h) (V) \g2) \z2it)J ' 

where (zi{t), Z2{t))£ H"^, and 

sup{t){\\zi{t)\\H^ + \\z2it)\\w^^) < e. 

t>0 



Remark 2.2. Certainly the global solution for p.Sp is also a global solution for <\2.1 
Our constructed solutions scatters in the iJ™ space at the optimal speed 1/t, i.e 



^v(t); ^ ' \f2j (V) V^ay H'" (t) 

Our method of proof also works for the 3D case. In that case we can obtain optimal 
decay 

Remark 2.3. To simplify the presentation, we do not make attempt to optimize the 
decay assumption (|2.9p on the final data (/i, gi, /2j 52)- More refined estimates can 
weaken the assumption p.9p significantly. However we shall not pursue this matter 
here. 

3. Reduction of the problem and the iteration scheme 
Define the linear solutions 

/ji(i)=cost(V)/i + 5^^||p5i, (3.1) 
, , sini(V) 

/l2(0=COSt(V)/2 + ^^52. (3.2) 

Let 
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For {u^^\v'-^'>) we have the system 

(□ + = V • ((v • V)(v(i) + /12) + uV(u(i) + /ii)) 

- 9t ((v • V)(u^i) +hi) + uW- (v(i) + /i2)) (3.3) 

(□ + l)v(i) - V((v V)(u(i) + /ii) + uV • (v(i) + /l2)) 

- dt ((v • V)(v(i) + /i2) + uV(u(i) + /ll)) 

-VA~^V(wv). (3.4) 



The RHS of p.3p - p.4p still contains quadratic source terms due to the self interac- 
tion of the linear solution {hi, /i2). To remove resonances created by these sources, 
we further define 

$1 - r '^^^^^^=fi^Gi{huh,){s)ds, (3.5) 



(V) 



$2= / 5i^l^i_f)^Z)G2(/^l,/^2)(s)ds. (3.6) 



00 



It is clear that 



,2 



< 



(□ + 1)$1 =Gi(/li,/l2), 
(□ + l)$2 = G2(/ll,/l2). 

Also we have the estimate 

ii<i>i(t)iiH. + miit)y\\m + w^-iimH^ + \\i>f2it)y\\H^ ^ . 

We then introduce one more change of variables and define 

u = - $1 = U - /li - $1, 

V = V*^^' - $2 = V - /12 - <i>2 . 

For the new variables (u, v), wc have the system 



(□ + l)u = V • (^(v • V)(v + $2) + uV{u + <Pi] 

- dt ((v • V)(u + $1) + uV • (v + $2)j 
+ V • (((v + $2) ■W)h2 + {u + $i)V/ii ' 

- dt (((v + $2) • V)/ii + (il + $i)V • /i2) , (3.7) 
(□ + l)v = V ((v ■V){u + $1) + uV • (v + $2) 

- dt ((v • V)(v + $2) + u\7{u + $1 
+ v((v + $2) • V/ii + {u + $i)V • /i2j 

- 9t (((w + $2) • V)/i2 + (u + $i)V/ii' 
-VA"iV(uv), (3.8) 
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and 

u = u + + hi, 

V = v + <^>2 + h2. (3.9) 

The system p.7p - p.9p in the variables {u, v) no longer contains resonant terms. 
We will construct a smooth global solution to this system. 

Let X e C^O^) be such that xit) = 1 for < ^ < 1 and x{t) = for t > 2. For 
integers n > 0, define 

Xnit) = X(2~"i)- 

It is clear that Xn is localized to the regime t <2^. 

Let = 0, v^"^ = 0. For n > 1 we inductively define 



X„(<)5t((v(""i) • V)(S("') +7.("-i'V • (v("))) 
V • ((v("-i) • V)$2 + m(""^'V$i) 
5t((v("-i) • V)$i + m("-i)V • $2) 

• (((v("-i) + $2) • V)/i2 + + $l)V/li) 

(((^("-1) + $2) • V)/ii + + $i)V • /la) , (3.10) 



- V 



and 



Here 



(□ + l)v(") = x„(t)v((v("-i) • V)(w(")) + V • v(")) 

- Xn(05t((v("-1) • V)(v(")) +«("-l)V(s("))) 

+ v((v("-i) • V)$i + u("-^)V • $2) 

- a* ((v("-i) • v)$2 + it*""^^ v$i) 

+ v((v("-i) + $2) • V/ii + (u("-i) + $i)V • /12) 

- 9t(((v("-i' + $2) • V)/i2 + (m^""'' + $i)V/ii) 
-VA-1V-(m("~iV("-i)). (3.11) 

Jn-l) =u("-l)+/ll+3>l, 
v("-l) = v("-l) + + $2- 



In Section Owe will study the system p.lOp ^ p.lip and its solvability. 
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4. SHARP PHASE ESTIMATES FOR $i, $2 

In this section we estimate (<I>i,$2)- By (|3?5 |) - (|3^ . (|3l l) -(|3?2 |) and (|2?2|) . we 
have 



$1 = ^Mt^^^iv) ((^^(^^ . ^^^^^^^ ^ h,is)yh^is)) 

- ds (/i2(s) • V/ii(s) + /ii(s)V • /i2(s)))ds, 

where we recall 

hiis) = coss(V)/i + -J^9i, 
h2{s) = coss(V)/2 + ^'"^^ 92- 

By Fourier transform, one can write down the explicit form of $1 in terms of the 
sum of several bilinear integrals involving the functions (/i, gi, /2, 52)- There are 
altogether sixteen such terms. To simplify the notations, we shall only write 

J 00 JR2 

where 

/(e)-(Cir^ + c2(0''^)/(0, 

and |ai| < 2, \a2\ < 2, ai, a2 are multi-indices, also ri — ±1, r2 = ±1, Ci, • • • ,6*4 
are constants, /, g can be any one choice from (/ii 5ii /21 52)- 
Let m > be any integer. We compute 

ll(t>$i(i)||H™+^ + ||(t>at$i(t)||^™+2 

where in the last step we made a change of variable 77 — > | + ry. 
The main result of this section is the following 

Theorem 4.1. For any integer m > 0, we have the estimate 

<csup(o™+^^(^ \d^frm+ J2 

^^^^ |q|<2 |q|<2 

1=1,2 1=1,2 



where C is an absolute constant, (/i, gi, /2, 32) are Fourier transforms of the final 
data (/i, gi, /2, 52) respectively. Same estimates also hold for $2- 

Theorem 14 . 1 1 will follow from the following 
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Proposition 4.2. For any integer m > 0, we have 

< sup (1^1 + \fi\r+'\\Hif,)\ + \id^h)iif,)\ + \{dih)i^,m (4.1) 

Note that it is enough to prove Proposition |4]2] for the case m = 0. We defer the 
proof of Proposition 14 . 2 1 to the end of this section. For now we assume Proposition 
14.21 holds and complete the 



Proof of Theorem \4-l\ It suffices to show the case m = 0. By Proposition 14.21 and 
the brief computation preceding Theorem 14. 1[ we only need to bound the quantity 

sup (III + \v\?\mM + \{d^mM + \{dih){im. 

where h{^,r]) ~ /(| + r])g{^ — rj), and f{S,), g{^) depends only on the final data 
(/i: 5ii /2, 32)- An elementary case by case analysis then yields the result. □ 



The rest of this section is devoted to the proof of Proposition 14.21 Wc begin by 
establishing some elementary lemmas. 



Lemma 4.3. For any a, b G M^, we have 

a + b a — b 



{a + b) {a- b) 



> 



(4.2) 



Proof. If a = or 5 = then (|4.2p holds trivially. Now assume both a and b are 
nonzero vectors. Denote the vector on the LHS of (|4.2p as X. Wc consider several 
cases. WOLOG wc may assume a • 6 > 0. 

Case 1: |a| > |6|. Then {a + b) ■ {a-b)> 0. Therefore we have 



{a + b)- {\a\ + \b\y 
Case 2: |a| < |6|. We decompose the vector a as 

a = ai + 02, where oi || 5, 02 -L b. 

Then we have 



y (a + 6) (a - b) [a + b) (a - b) 

Consider two subcases. 

Subcase 2a: |a2| > ||a|. Clearly by (|4.3p . we get 

lvl>l I 1 ■> W\ 
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Subcase 2b: |ai| > ^\a\. Since a ■ b > 0, we can write ai = Ab, with ^ ■ < 
C By (1131) and the fact ai 11 b, we obtain 

(A + 1)6 (A -1)6 



1^1 



> 



(a2 + (A + l)6) (a2 + (A-l)6) 



> 



(A+l)2|6|2 



(A-l)2|&p 



1 + |a2P + (A + 1)2|6|2 y 1 + |a2|2 + (A - 1)2|6|2 



(l + |a2p)-4A|fc|" 

> (l + |a2|^ + (A+l)^|b|^)-(l + |a2H + (l-A)^|b|^) 



> 



(A+l)"|bP , / (A-l)^lfcP 

1 + |Q2P + (A+1)|6|^ ^ V l + |a2|^ + (A-l)|&r-' 

(a-b) 



> 



(a + 5)2 . (a -6)2 
\a\ 



{H + \b\)^- 

The lemma is proved. 

Lemma 4.4. For any A, B eR, we have 



v/l + A2 + v/l + B2 - ^1 + (A + 5)2 > 



Proof. This is almost trivial. We compute 



1 



V1 + A2+V1+S2 



□ 



(4.4) 



LHs of Ea = 



2 + A^ + B^ + 2\/l + A2 • V1 + S2 - (1 + ^2 + ^2 + 2AB) 



x/TTA^ + VTTB^ +^1 + {A + By 



> 



1 



V1+A2 + Vl + 52 



□ 



Now define 



«i(c,^) = (e)- 

a2(e,r]) = (e) + 
a3(e,r])-(e)- 
a4(^,?7) = {0 + 



V) 
V) 
V) 
V) 



v), 
v), 
v), 
v)- 



Lemma 4.5. For any ^, 77 G M2, we /laue 
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< 



\v\ 



Vl< i < 4. 



1,2 



If i = 3, 4, then 



\d^M^,v)\<\^l Vl<fc<4, 



V|V,a.(C,r?)|^ 
|V,a.(e,77)|^ 



< 



< 



(4.5) 

(4.6) 
(4.7) 

(4.8) 
(4.9) 

(4.10) 



Proof. The lower bound in ()4.5p follows from Lemma 14.41 The upper bound in 
(|4.5p is trivial. If i = 1,2, then ()4.6p follows from Lemma [4.31 with a = ry, 5 = |. 
Similarly the case i = 3, 4 follows from Lemma l473l with a = |, 5 = 77. Next observe 
that \d';{{x))\ < {x)-<-''-^'> for any fc > 1. Therefore ^l^i holds. Next we show 
(j4.8p . This is a matter of direct calculation. Define a{x) = (a;). Then for i = 3,4, 
1 < A: < 4, 

\d^a.{^, v)\ = ma)iv + |) - (a^)(,7 - f )l 

<lel, 

where in the last inequality we have used again the fact that \d^{{x))\ < {x)^^''^^\ 
This finishes the proof of (|4.8p . Finally we show (|4.9p and (|4.10p . By direct 
differentiation, disgarding inessential coefficients and keeping only the main terms, 
we have 



( V„a,(^,7;) \ {d^a^{^,ri)Ydla,{(„ri) 



d. 



" VlVa.(e,'y)P 

/ Vr,a»(C,??) 

niva.(e,r/)p 



|V,a.(^,r;)|4 

g^.(C,^) ■ {^^a^{^,l^)y + (5,^a.(C,^))^ ■ d,Mi,ij) 
|V^ai(C,7y)|4 
(9,a.(e,7?))3 .(9X^,7?))^ 
|V,a.(^,7;)|6 



Therefore by (|46|) and (|48|) . we obtain 



d„ 



"V|V,a.(C,r7)P 



|V,a.(e,?7)P 
^ 1^1 ' 



GLOBAL SOLUTION FOR EULER POISSON 



17 



and 



< 



m + iv)' , ,.,2 m + 



1^1 



1^1= 



□ 



Lemma 4.6. Let aj(^, r/) be as defined as above. Define 



F^is,0 



„is(a.(C,ri)-ai(C,0)) 



h{ri)dr]. 



Then we have 



{■V 



\dsF.,{s,0\ 



\a\<2 



< 



|a|<2 



l';l<2 



|a|<2 



(E i\v\V'\d''Hv)\dv+J2 s^p|5"M^)l) 



|a|<2 



|r,|<2 



Proof. Wc first estimate 9^^";. Without loss of generafity, we also assume s > 1. 

We first consider dsFi, dsF2, where we shall use the following bound for a^, 
1=1,2, 



fivo-i > H 



(4.11) 



\d''ai\<l, forfc>l, 
[|a.(C,r?)-a,(e,0)| <|77|2. 

Since for the phase function a,;(^, 77) — ai{^, 0), the stationary phase point occurs at 
r] = 0, we introduce the cutoff function separate the near-origin regime and others. 
Let ipiri) be a smooth cutoff function such that 

'1, hl<l, 
0, |77|>2 



(4.12) 



We split dsFi, i ~ 1,2 as follows 

dsF{s,0= [ (a,(^,ry)- a, (C,0))e^^(''- («^o»/i(r;)dr7 

i 

(a,(?, v) - a.(C, 0))e^^(°-(«-'')-'^'(«^"»/i(r;)VX?7)rfr? (4.13) 

+ I (a.(C,ry)-a,(e,0))e'^('^'(«''')-'^-(«'°»;i(77)(l-V(?/))rf?? (4.14) 



For (|4.14p . we do integration by parts twice. This gives 



rV • 



■{a,{^,,j)^a,{^,0))h{,j){l-iir^)) 



drj. 
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Using the product rule and the estimates (|4.11l) . we obtain 

irai < 4 / (1 + \v\-'m\ + E \9"hiv)\)drj 

' •^I''I>1 |„|<2 
' •^I''I>1 |a|<2 



< 



\a\<2 



Wc now turn to estimating (|4.13p . To this end, we further write it into the foUowing 



For (|4.15p . we simply use (|4.1ip to get 



mB < / \v\\h{v)mv)\dv 

"'|?)|<2s-io 



|/i(77)|d77. (4.17) 



For (|4.16p , we use integration by parts twice and (|4.1ip to get 

< 4 1 / |a. (e, ^) - a. (^ , 0) \\Va,r'\h{^)Umi- i:{vs'°)\dv 



drj 



\a^{Ul) - a,{L 0)||Var'( E I^^^^WDdl " + s'°|V^(r,si")|)dr, 

l'?l<2 |„|<1 

|a,(^,'?)-a.(C,0)||Var'(E ■^"'^WD 

l'?l<2 |„|<2 

X (|1 - 1^(775^°)! + si°|VV(r?s'°)| + s'°|Ai^(77si°)|)dr; 
|Va,||Va,r3|/i(,y)|d^ 

hl<2 

+ / |Va,||Var'( E l^^'^WDdl - + s'°|V^(7?si°)|)d7y. 
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\mB\<\l [ \v\-'{\v\ + m''\h{v)\dv 



< 



|q|<1 

(H + ISI>'(E \9''hiv)\)dv 

|a|<2 

+ I^I)'(E \d''Hv)\>''drj\ 



M-^-'" \a\<l 



l"!---^" |a|<2 



^' |a|<2l''l<2 



Collecting the estimates above we obtain 



\dsF,{s,0\<^m'' TA (\v\V'\d"h(v)\drj+ snpld'^hm. (4.18) 



52-\i^i/ ^ 

|a|<2 



|r,|<2 



for i = 1,2. 

Now it remains to obtain the estimates for i^i, i = 3, 4. We will use the following 
estimates for Ui, 



■|va.|>WnW- 

\d''a,\<l, k>l, (4.19) 
|a^(e,^?)-«.(^,0)| < |77|. 
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We will estimate dsFi simply by integration by parts twice. We have 



\VaA 



Va,\~^h{TjMdi^ 

+ / |Var' E |a"/iW|(h| + |Va,|)d77 
l«l<i 

|a|<2 ^ 
^ l"l<2 



|a|<2' 

This finishes the estimates for dsFi, i = 3,4. Finally we make a remark about 
the decay of Fi. Indeed, the decay of Fi can be obtained by slightly modifying 
the above computation, namely, replacing the integration by parts twice with IBP 
once. We omit the details. The proof of this lemma is completed. 

□ 

As a consequence of this lemma, we have 
Lemma 4.7. Assume the function h = h{^,T]) is such that 



Then for i = 1, 2, 3, 4 and any ^ 7^ 0, we have 



Proof. Define 



K2 iai{^,r]) 



h{^,r])dr]. 



(4.20) 



We write 

/•OO 



1 



1 



1 



m,(e,0) 



za,(C,0) 
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Now due to the decay of dsFi, the integral in the above is absolutely convergent. 
Therefore we have 



e"''^<^i'''^hi(,ri)dr]ds 



= lim 
= lim 
= lim 

£^0 + 

lim 



ia{£,, 0) — e 



h(^, ri)dr]ds 



h{^,r])dri 



Lemma 4.8. Let Q < 6 < . Then for i = 1,2, 



j 



drj 



□ 



(4.21) 



Proof. This is a stationary phase calculation. The main point is to get the explicit 
dependence on the parameter ^. WOLOG we can assume t > 1. Define A = 
By rotation invariance, we may assume ^ = 2Aei = (2A,0). Let rj = Xfj. Then 

(^ + f) + ('?-f> 



= Vl + A2(|?7P + 2771 + 1) + + A2(|7?|2 - 2771 + 1). 

Here ?7i is the first component of the vector ij. Passing to radial coordinates, we 
obtain 

\/l+4A 



LHS of Km < 



2tt 



(,ita(pfi)x^pdpde 



where 

g{p,e) = v/1 + A2(p2 + 2pcos6l + + + X^{p^ - 2pcose' + 1). 
Now set r = y/p. Then we get 

r.27r 



LHS of iKm < 



(4.22) 



Here 

~g{r, 9) = y^l + A2(7' + 2\/7cos6' + l) + 1 + \'^{r - 2^^0080 + 1) 
Observe that for 7' < jq^p-, 

1 < A+, A_ < A + 1. 
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On the other hand 

d ^, ^2 A+ + A^+2^A+-A^-4X^cos-. 
— q[r,0) = \ ■ , ; ; 

2 2 + 2A2(r + 1) + 27(1 + A2(r + 1))2 - 4A4r cos2 6 - 4A2 cos^ i 



A2. 



For r < jq^p- , we have 



4AVcos2 6' < 2A2 < 2A2(r + 1). 



This gives 



2 + 2A2(r + 1) + 27(1 + A2(r + 1))2 -4A4rcos2 6' - 4A2 cos2 
>2 + 2A2 + 2A2(r + l)-4A2 
>2. 



Therefore we get for < r < j^^, 



> 



A2 



(4.23) 



(1 + A)t 

Next we need to bound the second order derivative -^g{r,9). We first rewrite 
g{r, 9) as 



ff(r,0) = y^(717+ 717)2 

= ^^2 + 2A2(r + 1) + 27(1 + A2(r + 1))2 - 4A4r 



cos2 9 



= : 7^1 + 27^2 
By direct differentiation, we have 
92.9 1 



- • (Si + 27B^)-^ • (--^2 = • {-^f + i?2 



It is not difficult to check that for < ?' < 



Bi+27b^>1, and B2 > 1, 



Qj.2 



dBi 



dr 



Also 



Therefore 



< A2, and 



92^2 



gj.2 



dr 



< A2 + A^ 



925 



Qj.2 



< (A2 + a4)2 + 

< A^ + A*. 



(4.24) 
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By ((423| and (|424|) . we obtain 



IBS-' 
1 + 4A^ 



< _ 



< 



1 ( {1 + X)i , (1 + A)3 
A2 

1 /(1 + A)i 



A2 



A4 

+ (1 + A)5 



(A^ + A«) 



1 + A2 



Therefore by dOU, we get 



1 



LHS of Km < J ■ (1 + A)^. 
The lemma is proved. 

Lemma 4.9. Let -ijj e C^'{R^) be such that 

^ ^ \0, |x| > 2. 
T/ien for any i ~ 1,2, we have 



l')l< 



1 



<--.(0'- sup (IMC,^)l + |v,Me,'7)l)- 

w lf,l<li# 

Proo/. If t < 1, then ((4:251) holds trivially. Now assume t > 1. Define 



Note that ipiO) = 1 and 



By Lemma 14.51 we have for i = 1,2 



1^1 



sup \G,i^,v)\<{0 ■ sup 



l'7l<' 



Since Gi(^,0) = 0, we get 



sup ^'^^P |V'^G,(e,77)| 

l^l<^ l'7l 



l'7l< 



< 



sup 



|V,/i(^,r?)| , |Me,r/)|-|V,a.(e,7?)| 



l„l<21i£l V |a«(C,'7)l 

,(0^ sup l'^^^'^)" 



1^1 



hi' 



25|£| |a,;(C,T?)| 



< 



(0! 



sup (|/i($,r;)| + |V,Me,r;)|). 



I';l< 



Tip" 



□ 



(4.25) 



(4.26) 



(4.27) 
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Using integration by parts together with ()4.26|) . (|4.27p . we then obtain or i = 1, 2 
LHS of 6221) = i| / G,i^,rj)^^^-fje^'-'(^'^^da{v) 



V 



1 r 

< 



^ ^ (0- sup \hi^,ij)\.M±Mlda{^) 



+ sup (|MC,^)l + |V,Me,^)l)- / <'^' + '^'^' dry 

<i.(0'- sup (IM^,^) + |v,Me,^)l). 

□ 



Lemma 4.10. For i = 1, 2, 



have 



1 



<-.(l + |log|e||)-(0-'- sup (||| + |7?|)^^-(|Ml,^)l + |V,Ml,'7)l). (4.28) 



is localized 



Proof. WOLOG wc may assume t > 1. Since the function 1 — V'( 
to the region {\r]\ > and vanishes on the boundary, a simple integration by 



parts gives us 



LHS of gSHD = -- / V 



(l-^(^))) e^'^^^^'^U^ 
(4.29) 



It only remains for us to bound Gi(^,ry). 
By direct differentiation, we get 

|G.(e,.)i<#4H^i-^-|i-^(-'^^^^'' 



|V,a.(C,r])P |a,(e,?7)| 
|V^a,(e,77)l ' |a.(^,7?)| 

|i-V(- 



|i-V'(^)l 



|a.(e,r?)P 
1 



6\^\ " 



|V,a,(C,r?)| |a,(e,?7)| <5|e| 



(4.30) 
(4.31) 
(4.32) 
(4.33) 
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By Lemma l4.5|, we have 



< m + \v\y 



|^|2 IM^,^)l-i|,|>^ 



Similarly we compute 



&<m + \v\)^-\hitv)\-\,yliii, 

6331) < — ^ • ^ • ■ 

Collecting all the estimates, we get 



sup (||| + |ry|)^.(0^.(|Me,'7)l + |V^Ml,^)l) 



■ sup {\i\ + \v\y'-i\hiifj)\ + \y,hiim- 

Plugging the last estimate into (|4.29p and integrating out 77, we obtain (|4.28p . □ 
Lemma 4.11. Let i = 3,4. Let t > I and \^\>j. Then 



^T^-JTr-lh- (l^l + l^l)''(IMl,'?)l + l(5,ft)(|,^)| + |(9|/i)(|,77)l). (4.34) 

Proof. By using integration by parts twice, we get 



V2 a^{^,v) 



^ / G'2(e,77)e'*'''(«''')dry. (4.35) 
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It only remains for us to bound G2{S.,r/). By direct differentaion and discarding 
inessential constants, we have 



< 



"VlV,a,(e,,7)P 



V,,aj(^, 77) 



1 



1 



|V„a.(C,?7)| 
1 

1 



" V|V,a,(C,ry)| = 



|V,a,(e,7y)P 
1 



|V^a.(C,??)| |V^a,(C,77)l 
For (|4.36p . we use Lemma [4.51 to bound it as 

(|4-36P < 1^^^ \h{i,m 



di 



(4.36) 
(4.37) 
(4.38) 
(4.39) 
(4.40) 



Similarly 



< 



< 



< 



'''V|V,a,:(e,ry)P 

/ Vr)a»(C,??) 

"'MV,a,(^,,y)P 

-I- lnl\8 

1 



1 



1 



1 



a»(^,'7)l |V,,aj(^,77)| 



n\r(.\h{lfi)\ + \id,h){iv)\). 



1 1 

The computation of (|4.38p - (|4.39p is similar to (|4.37p . and we have 
Finally for (|4.40p . wc note that by Lemma 



Htv) 



< 



\id'h)i^,v)\ 



\i^nhm,v)\■\^^a^{^,v)\ 



< 



(lei + + \id,h)i^,v)\ + l(a»(^,^)|) 



Therefore 

< 



1 1 1 



i\h{tv)\ + \id,h){^,v)\ + \id',h)i^,v)\) 



i^p (0^ ivr 



■ sup ( 



'mifi)\ + \idr,h)iifi)\ + \idlm,m- 
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Collecting all the estimates, we obtain 

Is I \5/ W) i,f^es.^ 

Plugging this last estimate into RHS of (|4.35p and integrating out ij, wc obtain 
(|4.34p . The lemma is proved. □ 



We are now ready to complete the 



Proof of Proposition \4-S\ Assume h = h{£^, if) is such that the RHS of (|4.ip is finite. 
By Lemma Wl\ we have for all 1 < i < 4, 



Consider first i = 1, 2. Let e C;f (R^) be such that 



ip{x) 



1, |2:|<1, 
0, |x|>2. 



Let < 5 < Then we decompose 

gitai(5,r)) 



(4.41) 



drj 



(4.42) 
(4.43) 



For (|4.4ip . we use Lemma [4. 101 and this gives us 



Jr2 ai(?.^) 



< 



sup m + \v\y' ■ + \vnHc,m -m + i logici) • (o^ 



< sup {\i\ + \fj\)^^.{\h{lf^)\ + \V^h{i,f,)\). 
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Similarly for ()4.42|) . we use Lemma l49l to compute 



/A\2m / i Htv) ^MO\ hiC,0) \ 



< 



< 



(0" sup im,fi)\ + \Vfihi^,m 



snp m + mv' ■ {m,fi)\ + iv,m?,^)I) • iko^'l? 



< sup (||| + b?|)"-(IMe,^)l + |v,Me>7)l). 



For (|4.43p . we apply Lemma Lemma [4.51 and obtain 



ryitai 



< 



(0 



9 Hto) 



< 



my°hic,o)\\L-i 



sup m + \ri\y^-miv)\ + \Wfihii,m. 



< 

Therefore we have shown that in the case i = 1,2, the contribution of is bounded 
by the RHS of (|4.ip . Next we consider i = 3,4. Ifi<l, then by Lemma |431 we 
have 



E2 aj(C,??) 



< 



< sup (||| + |ry|>i2|Me,^)l 



which is acceptable. If t > 1 and |^| < i, then we also have 



< 



(Ifl + H>IMC,^)M^ 



i^(l?l<i) 



< sup {\i\ + \v\)''\h{lv)\. 
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It remains for us to consider the case t > 1 and |^| > j. By Lemma 14.111 we 
compute 



< sup m + \v\rH\Hiv)\ + \id,h)ii,^)\ + mh)i^,v)\) 



1 1 1 



< sup {m + \v\r\\h{lv)\ + \{d^h)ii,f,)\ + mh)ilf^)\). 



This finishes the case i 



3,4. 



□ 



5. Solvability of the iteration system p.iop - p.lip 



Since the RHS of p.lOp - p.lip generally contains second order derivatives of 
v^"^), the solvability of the system is not immediately a trivial matter. How- 
ever, due to the time cutoff function Xn{t) which vanishes for t > 2"+^, these 
troublesome terms involving (m^^^v^")) also drops out for t > 2"+^. Our strategy 
then is to solve first p.lOp - p.lip for t > 2"+^. In this regime we can make use 
of the explicit representation of the solution in terms of the RHS of (|3.10p - (|3.1ip 
which only involve ($i,$2j ). We then solve backwards the 

same system for t < 2"+^ with initial data (u(")(2"+i), v(")(2"+i)). In the latter 
case, the time interval is of finite length and the solvability is rather obvious since 
it is a linear equation for (w^"', v^")). 

More precisely, we inductively assume that 



("-i)(i))'|k.. + |lv("-i)(i)|U™ 

+\\{^^"-'Ht)y\\H^<-^. 



(5.1) 



Then for t > 2"+\ 



we get 




+ V • (((v("-i) + $2) ■y)h2 + («("-!) + $i)v/ii) 

" 9t(((«^""'' + *2) • V)/ii + (i2("-i) + $i)V • 



(5.2) 
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and 

(□ + l)v(") = v((v("-i) • V)$i + it("-i)V • 

- at((v("^i) • v)$2 + u^"~^'v$i^ 

+ v((v("-i) + $2) • V/ii + + $i)V • /i2 

- 9t(((v("-i) + $2) • V)/i2 + (w^""'' + $i)V/ii^ 
-VA-iV-(w("-iV("-i)). (5.3) 



Observe that the RHS of (|5.2p docs not contain u'^"^, therefore we have the 
exphcit representation 

= r ^^^^^^7^^Gl($l,$2,/^l,/^2,s("-^^v("-l))(s)ds, (5.4) 

for t > 2"+\ 

where 

Gi($i,$2,/^l,/l2,^l("-'\v("-l') 

= V • ((v(""i) • V)$2 + w("~^'V$i) (5.5) 
_ a* ((v("^i) • V)$i + V • $2) (5.6) 

+ V • (((v("-i) + $2) ■V)h2 + (m*"^^' + $i)V/ii) (5.7) 
- at(((5("-i) + $2) • V)/ii + (s("-i) + $i)V • /12) . (5.8) 
Similarly for v^") (t) we have 

v(")(0 = f G2($i, $2, /.2, ^("-^\v("-i))(.)d., (5.9) 

for t > 2"+\ 

where 

G2($1,$2>1,/12,S("-'\V(""1)) 

= V (^(v("-i) • V)$i + V • $2) 

- dt ((v("-i' • V)$2 + u'""^' V$i' 
+ v((v("-i) + $2) • V/ii + (w("-i) + $i)V • h2 

- a* (((v("-i) + $2) • v)/i2 + (s("-i) + $i)v/ii) 

^ vA-iy- (w("-iV("-i)). (5.10) 

For the function ($1, $2), we have 

||$i(s)||h™+2 + ||$2(s)||h-+^ + ||($i(s))'||j^,„+i 

+ ||('f2(.s))'||ff,. + l < V5>0. (5.11) 
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Also for (/ii,/i2), we have the estnuate 



\\hlis)\\L'^ + \\d,h,{s)\\L^ + \\D"'+^h,i.s)\\L 



dsD"'+^h^{s)\\L^ <-^, Vs>0. (5.12) 



e 

\\h2{s)\\Lr + \\dsh2{s)\\L^ + \\U^+^h2{s)\\L^ 

+ \\dsD'^+^h2{s)\\L^<-^, Vs>0. (5.13) 
(s) 

Recall that 

Therefore by dSU), ((5lT|) - (|5A3| . we have 

||V • ((v(""i'(s) • V)$2(s) + ^.("-i)(s)V$i(s)) 
<||V • ((v("-i'(.?) ■ V)$2(s))||ff™ + ||V • (s("-i)(s)V$i(s))|U™ 

+ ||V • ((/l2(s) • V)$2(s)) llff™ + ||V • (/li(s)V$i(s))||ff™ 



+ ||V • (^($2(s) • V)$2(s)j IIh". + ||V • ($i(s)V$i(s))||a™ 

<ll(v("-'H^))'lk". • ll*2(=s)i|H,„+. + ||(ft("-l)(s))'||«„. . ||$i(.)|U„ 
+ {\\h2{s)\\L^ + P"+1/i2(.s)||ls^) • ||$2(.s)||h™ + 2 

+ {\\hi{s)\\L^ + W^+'h,{s)\\L^) ■ ||$i(,s)||^™+2 

+ ll<I'2(.s)||?,„ + . + ||<i>l(s)||^,„ + . 



The remaining terms (|5.6p - (|5.8p can be estimated in a similar fashion. Collecting 
the estimates, we get 

||Gi(<i>l,<i>2,/ll,/l2,s("-'\v("-l))(s)iU.. < -4^. 

By (|5.4p . we then compute 

/C30 
|lGi($i, $2, /ll, /i2, 12("-1\ v("-l))(s)||ff^ds 

.oo g2 

Vt>2"+i. 



Similarly 



||s(")(t)||H™ + \\d.,M''\t)\\H^ < ^, V< > 2"+\ 
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Therefore 



£2 



+ < Vt > T+\ (5.14) 

By using similar estimates, we obtain from (j5.9p 

l|v(")(t)||H". + ||(v(")(t))'||H". V<>2"+\ (5.15) 



The estimates (|5.14P " ()5.15p are consistent with the inductive assumption (|5.ip 
on the time interval [2"+^,oo). However we still have to solve p.lOp - p.lip for 
< t < 2"+^. Since we have obtained the solution on the time interval [2"+^, oo), 
we can use (m*^"^ (2"+^), v(") (2"+^)) as initial data which is explicitly given by ()5.4p - 
dSll), i.e. 

^("^2"+^) = f ^)(^) G,i<l>,, h,, h,, v("-^))(.)d., 

-Jn)(^n+l, r Sin(2"+1 -S)(V) - („_i) („_,)^ 



-t„j(2"+^)= / ^ " ^ G2($l,$2,fel,fe2,^I^"-^^v^"-^0(s)dg. 

Using this information, wc then solve backwards (|3T0)) - (|3lT1) until i = 0. This 
way we have obtained the solution (M'^")(i), v("^(t)) well defined for all t > Q. 
Furthermore since [0,2"+^] is a finite time interval and ('u*^"^(i), v'-"-'(i)) enjoys the 
decay estimate (|5.14p - (|5.15p . we have 

|ls(")(i)l|H"> + ||vW(t)|j^,„ + ||(s(")(0)'|U"> 

ll(v("'(i))'llH."<|j^, Vt>0, (5.16) 

where Kn is a constant depending on the iteration step n. Because of the depen- 
dence on n, the estimate ()5.16p is certainly not good for us since we have to justify 
the inductive assumption ()5.ip for step n. Nevertheless it is an important a priori 
estimate for our argument later. In the next section wc shall upgrade the estimate 
(|5.16p to the true decay estimate (|5.ip . 



6. A PRIORI ENERGY ESTIMATES OF THE ITERATION SYSTEM p.lOp p.lip 

In this section wc carry out basic energy estimates for the iteration system (j3.10p - 
p.lip . Our standing inductive assumption is 

+ W'^''-'\t))'\W-<-^- (6.1) 

For the functions ($i,$2), recall the estimate 

||<i>i(t)|U,„+2 + ||$2(i)||H™+^ + ||($i(t))'|U™+i 

+ ||($2(t))'||ff.>^^i Vt>0. (6.2) 
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Also for (/ii, /i2), we shall use the estimate 



WhiimL^ + \\dth,it)\\L^ + \\D"'+'h,it)\\L^ 

+ \\dtD''-+'hiit)\\L^ <-^, Vt>0. (6.3) 

ii/i2Wiil- + \\dth2mL^ + \\n"'+^h2mL^ 

+ \\dtD"^+^h2mLr<4^, Vt>0. (6.4) 

{t) 



Our task is to justify the estimate (|6.ip for (m("',v(")). For the convenience of 
notations, we introduce 



+ ll(v("'(t))'llH™, Vn>0. 



Take derivative on both sides of p.lOp and multiply by the factor dtD^u^"^^ 
Integrating by parts and summing over |a| < to, we obtain 



I Q:| <m 



'-Xn{t) y / D^h- ((v("-i) • V)(v(")) + u(""i)Vu("))) dtD^u^'^^dx (6.5) 

-xn{t) [ £'"atf(v("-^) •v)(M(")) + u("-^)v-(v(")))ati:i"iz(")dx 

(6.6) 

+ y / i:>"V- ((v("-i) •V)$2+u^""^^V$i)9ti:»"u(")rfa; (6.7) 

\a\<m 

-V / i:i"at((v("-i) • v)$i + ^("-1^ • ^alatD^u^")^^ (6.8) 

\a\<m 

+ y / i:>"v- (((v("-i) + $2)-v)/i2 + (u^""^' + $i)v/ii)at£i"u(")dx 

|a|<m 

(6.9) 

y y" £l"at(((v("-i)+$2)-V)/ii + (u("-i)+$i)V-/i2)ati:i"u(")da;. 



|a| <m 



(6.10) 
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We first deal with ()6.5p . By using integration by parts and Leibniz in time, we 
write 

5l = Xn{t) V / L'^V- ((v("-i) •V)(v(")) + u("-i)Vu("))ati:'"u("'rfa; 

\a\<m — l 

(6.11) 

+ |(x«W E J^yr'^D'^d^^d^M'^^D'^u^-^dx) (6.12) 



H-- 

l<i,j<2 



l<jj<2 

I Q I —m 
l<i.j<2 



(6.13) 
(6.14) 
(6.15) 



\a\—m 
l<i,j<2 



By using (|6.ip . (|6.2p - (|6.4p and Sobolev embedding, we compute the absolute 
value of (|6.1ip as follows 



| (IOT|) | < E E / l^^v("-l'| • |i^"+2-/3^(„)| . 

|Q|<m-l |0|<|q| + 1"^*^ 

+ / \D^u^''-^'^\-\D°'+^-^u^''^-\dtD"iS'''^\dx 

^ E E IP''/i2(t)||Ls=-ll^"'''"^vHWIlLrll5*^""'"'WIU^ 

|a|<m-l |^i|<|Q| + l 

+ E E ii^^v("-i)(t)iu^. .|p"+2-^v(")(t)iu..|i5*i5"^2(")(t)iu. 

|Q|<m-l |^J|<|q| + 1 
l/9|<^ 



+ E E ll^^v("-iHt)IUr ll^"^'"''v(")(i)iu^. .||a,i?"K(")(i)IU^ 

|Q|<m-l |/3|<|q| + 1 

< 

The term (pT^ will be absorbed into the LHS of ([13]). To bound ((O^ . note 
that 

IxUOl < 2-"|x'(2-"t)| < 1, Vn>0, t>0. 
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By using integration by parts, ()6.iP " (|6.4p and Sobolev embedding, we obtain 

I Q I —771 

< {\\h2mLr + \\Dh2mLr + ll*2(t)||ff." + \\^^''-'\t)\\H^)Ar.{tf 

Proceeding in a similar manner, we also obtain 

For (j6.15p we shall integrate by parts twice and this gives us 

, I JR^ 



\a\—m 
l<i,j<2 



I a I —771 
l<«,i<2 



(6.16) 
(6.17) 
(6.18) 



I a I —771 
l<i,j<2 



We shall keep the term (j6.16p since it will cancel with the corresponding term 
obtained in the calculation of v*^") later. The term (j6.17p can be bounded after one 
more integration by parts in dx^ and this will give us 

\m\<-^A„{tr. 

The term (|6.18p can be bounded in a similar way as (|6.13p and we get 

\mB\<-^A„itr- 

Collecting all the estimates, we obtain 

\a\=ra ^ 
l<i,j<2 



at) E / 



H-- 

l<i,j<2 



where 



\E,\<4-A„{t)\ 



(t)' 
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Next we consider ()6.6p and write it as 

(6.19) 
(6.20) 



I a I < rn — 1 * 



tn(i) y / i:'"9tf(v("-i) • v)s("')ati:'"u(")(ix 

|a|— m 



(6.21) 



The term (|6.19p can be easily bounded by 



For (|6.20p . by using integration by parts, (|6.ip - (|6.4p . Sobolev embedding and 
proceeding in the similar way as in the estimate of (|6.5p . we obtain 



[ ((Om i < y 1/ i?"((atv("-i) • V)?2("))9tZ?"?2(")d2; 

<y y / |9ti:»"-'3v("-i)| • |vi:»'^iz(")| • |9ti:i"u(")|d2: 

|a|=m|/}|<|a|-^»' 

y I / (v("-i) • V)^"^*^'") • dtD"u^''Ux 



\a\—m 



For (|6.2ip . we use Leibniz in time and write it as 

(lOTj) ^ - y / i:)"(atw("-iV •v("))9ti:)"i2(")dx (6.22) 

\a\—rn 



(6.23) 



\a\—rn 

We can easily bound (|6.22p as 



|Q|=m|/3|<|Q| 
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On the other hand for ()6.23|) . we further decompose it into 

|a|— m 

|a|=mO<|/3|<|Q| -^"^ 

= y / u("-i)5tL'"v(") • ydtD'^u^'^Ux (6.24) 

|q|— m 

+ V / (vu("-i) •ati:i"v("')at£i"M(")dx (6.25) 

|ck|— m 

- y I Df^u^"~^')D"-^dtV ■v^"^dtD°'u^''Ux. (6.26) 



We shall keep the term (|6.24p since it will cancel with a corresponding term 
arising from the calculation of v'"-' later. The terms (|6.25p . (|6.26p can be bounded 
in the same way as the calculation of (|6.5p and this gives us 



\mB\ + \m^\<^Mt)'- 



it) 

Collecting all the estimates, we obtain 



+ E2 



where 



\E,\<-Mtr. 



Wc still have to bound the terms (|6.7p - (|6.10p . The estimates of these terms can 
be done in a similar manner as that of ()6.5p . Wc omit the tedious computations 
and write the final result as 



\m\ + \m\ + \m\ + i»i < ^Mt)- 
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Putting together all the estimates, we have the expression 



■{Y, \\D"u^-\t)\\l.+\\(D'^u^-\t))'\\l: 



\o(\<m 



I Q I — m 



l<i,j<2 



where 



F1+F2, (6.27) 



\F2{t)\ < 77T2^nW- 



Next we consider the bounds on v^"^ by using p.lip . Proceeding in a similar 
fashion as in the computation of u^"'\ we obtain 

^l(E ll^"v(")(t)||i.+||(W")(t))'||i; 
\a\<.7n 

=x 



where 



|a|— m 
l<i,i<2 

+ + F4, (6.28) 
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Denote 



|Q|<r 



J2 ll-D"v(")(t)||2, + j| ('i:)"v(")(t) 



2 



2 



|a|<m 



\a\—m 

By using (|6.ip - (|6.4p . we have 

|(imi + |(|S3ni)| + |(IS3Il)| 

<e.i( ^ ||^"s(")(i)||i 



(t) J2 J 9,^vj""'^I?"9,^vf^D"u(")dx (6.29) 

|a|— m ^ 
l<i,j<2 

(t) V / v^""'^D°5,,vf^9,^D"M(")(ix (6.30) 

|a|— m 
l<ij<2 

(i) y / w("-i)|i?"V- v(")|'dx. (6.31) 



I Q I < m 



.+11 (W")(i))'| 

5^ pV")(t)||i.+||(W")(t))'| 



2 

Li 



|Q|<r 



Therefore for e sufficiently small (depending only on an absolute constant), we 
have 



3.32) 



where C is some absolute constant. Adding together (j6.27p and (|6.28p and using 
we have 

^2 

(6.33) 



|i?«w|<^i?„(t) + |„ t>o, 



where we have canceled a factor of Bn{t) on both sides of the inequality. 
By (|5.16p and (|6.32p . we have the rough estimate 

\Bn{t)\<^^ yt>o, 

where Kn is a constant depending on the iteration number n. 
We now need the following 



(6.34) 



Lemma 6.1 (Infinite time Gronwall lemma). Let B(t) he a nonnegative smooth 
function on [l,oo) satisfying 



K 



B{t)<—, Wt>l. 



(6.35) 
(6.36) 
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Here < e < 1, (5 > 0, A' > are constants. Then 



B{t)< 



1 



1 - e t 



-, Vi>l. 



3.37) 



Remark 6.2. The crucial point here is that the final bound on B{t) is independent 
of the constant K . 

Proof. By Fundamental Theorem of Calculus and plugging the estimate ()6.36|) into 
(j6.35p . we obtain 



B{t) < 



dT 



B{t) 



dT 



< 



+ ^]dT 



< 



T T 

Ke + S 



t 



We inductively assume j > 0, 

B{t)< 



AV+i + (5(1 



Then for j + 1, using (|6.35|) we have 



dt 



Bit) 



< 



t^ 



e{Ke 



5{l 



Using again the fundamental theorem of calculus we conclude that the inductive 
estimate holds for all j . Taking j — > oo we obtain (|6.37p . □ 

By Lemma WA\ and (|6.33l) . (j6.34p . we obtain for e sufficiently small 



Bn{t) 



< 



3.38) 



The extra power of e can be used to kill the underlying implied constants. Therefore 
the inductive assumption (|6.ip also holds for step n with the same underlying 
constant. Hence we obtain the following: there exists a constant C > 0, such that 
for all n > 0, 



+ll(s(")(t))'llH™ + ||v(")(t)|U™ 
+II(v(")W)'||h"><c^, yt>o. 



(6.39) 



Remark 6.3. Since the bound ()6.38p is slightly stronger than the inductive bound 
(|6.ip . one may wonder whether the bound (|6.ip is optimal in terms of the power of 
e. Indeed by revisiting the inequality (j6.33p and Lemma l6.1i one can easily derive 
the bound 
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from which it follows that 

|ls(")(t)|U.. + ||(i2(")(t))'||H™+||v(")(t)||H™ 
+ m^"Ht)y\\H^^<C^, Vi>0, n>0. 



7. Contraction estimates and existence of the smooth solution 

In this section wo shall show that the sequence (u^"-', v'"^) forms a Cauchy se- 
quence in some suitable metric. Define 



[/(") = y(") _ 



From ([Q])) . (i3lT|) . we get 



- Xn{t) ■ dt ((v("-2) . V)^(n) ^ y(«-2)y . -^/(«)^ 

+ V • • V)$2 + C/'""^V$i) 

- • v)$i + [/("-i)v • $2) 

•V)/n + ^'V -/la), (7.1) 



and 
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+ Xn(<)v((v("-2) . V)([/(n)) + y("-2)v . 

- ((1^^"-'^ • V)(v(")) + [/("-1)V(S("))) 

- X«(<)<9t ((v("~2) . v)(f(")) + u("-2)v([/("))) 
-(XnW-Xn-l(^)M((v("-2)•V)(v("-l)) + ^i("-2)V(w("-l))) 
+ v((T/("-i) • V)$i + [/("-i)V • $2) 

- 9t((F("^i) • V)$2 + [/^"-i)V$i) 

- dt{(y^"-^^ ■ v)/i2 + ?7("~iV/ii) 

- VA-ly • (7 2) 



For the convenience of notations, we define 



+ ^>0. (7.3) 



We shall derive some recursive inequalities for Wn{t) from which it will follow 
that the sequence (u("\v(")) is a contraction. 
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Multiplying both sides of ()7.ip by dtU'^^^ and integrating by parts, we obtain 

Xnit)W ■ f • V)v(") + U^"-^^\7u^"A ■ dtU^"'^dx (7.4) 

2 V / 

+ Xn{t) I V • ((v("-2) . y^yin) ^ ^(n-2)^^(n)\ . Q^Jjin)^^ (7 5) 

+ (Xn(i) - Xn-lit)) f V • ((v("-2) . V)v("-1) + ^("-2) y - („- 1) \ 5^ rf^ 

(7.6) 

- Xnit) J dt • V)it(") + f/("-i) V • v(")) atJ7(")dx (7.7) 

- Xn{t) I dt f (v("-2) . V)t/(") + li("-2)y . ^/(") [/(") rf^ (J.g) 



(7.9) 



(7.10) 



+ / V • ((V^("-i) • V)$2 + C/("~''V$i)atC/(")da; 

dt f (^("^1) • V)$i + [/("-i)v • $2)5tf/(")da; (7.11) 

V • ((^("-1) . V)/i2 + t/("^i)V/ii)atC/(")dx (7.12) 

5Y(V^(n-i) .V)/ti + t/("-i)V-/i2)<9tt/'"'da;, (7.13) 

By using (|09)) . (I01)-(l63) and ([LSI), we have 

I (izai) I + 1 (Ell) I + 1 Gini) I + 1 Gill) I + 1 

Also we have 

<\Xn{t)-Xn^l{t)\--^-Wr,{t). 



By using integration by parts, we write (|7.5p as 

(1731) = -x„(t) / u("-2)V[/(") • dt^U^^^Ux (7.14) 

^(*) E / v("-2)a,,F/")5ta,^.f/(")dx. (7.15) 



Xn 

l<i,j<2 ' 
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For ()7.14|) we use Leibniz in time and this gives 

(Ell =-^(ix«W^^"'"-'^|VC/("'pdx) (7.16) 

+ lxn{t) I I V[/("fdx (7.17) 

+ \xn{t) I 5tii(""2)|V[/(")|2dx. (7.18) 

The first term (fTTel) will be absorbed into the LHS of (fTi)) . By using ((09l) . 
(j6.2p - (j6.4p . the other two terms ()7.17p . (j7.18p can be easily bounded as 

\m\\ + \^\<j^^n{tf. 

Next we consider the term (|7.15p . By using repeatedly integration by parts, we 
have 



l<ij<2 



it) E /^9..v("-2'9,^y/")9tC/(")dx (7.19) 

-X«(t) E / 9,^^,.vf-'V;"^atC/(")da; (7.20) 

+ Xn(t) E /, v!""'^9.,y/"^9.,C/(")rfa:) (7.21) 

+ x;(0 E / a..vl"-^V/")9..C/(")dx (7.22) 

l<i.j<2 -^^^ 

+ xUO E / v("-^)a.^.F/")a..[/(")dx (7.23) 

+ X«(t) E f ^d,^dM''~^^vl'''>d,,U^"Ux (7.24) 

+ Xnit) E / atv("-')a,^^/")a,,C/(")dx (7.25) 

l<ij<2 "^^^ 

-Xn(<) E / 9t^/"^5,^(5.,C/^"^-vf-'))dx (7.26) 



l<?,j<2 * 

+ Xnit) E / aiF/"^5,,t/(")9,^.vf-')dx (7.27) 

We shall absorb (|7.2H) into the LHS of (j6.39p . We will make no changes to the 
term (|7.26p since it will cancel with a corresponding term in the calculation of V'-^^K 
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For the rest of the terms, we use again (|6.39p . ()6.2P " ()6.4p to bound them as 

|(i7i^l + 1(17:^1 + 1(17:^1 + 1(17:^1 + 1(17:^ 

This ends the estimate of (|7.15p and ()7.5p . By a similar but tedious calculation, 
we also have 



xn{t) [ dtU'^-^^yu^^A ■dM^'^dx 



where 



\Ei\<^Wn{ty. 
Collecting all the estimates, we obtain 

^{l\\u^"Hmk + l\\{u^''Hmk 
d /I 

~ di 



l<ij<2 -^^^ 

.lo2 V / 



E2, (7.28) 



where 



1^2 1 < -Kw;„(i)' +u.„(i) • u;„_i(t)) 



+ IXnW - Xn-1 • J-^Wn{t). 

Next we deal with the estimate of y^"). This time we need to take the inner 
product of both sides of (|7.2p with dtV'^"'^ integrate by parts. Proceeding in a similar 
manner as in the calculation of [/'■"•' and after a long and tedious calculation, we 
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obtain 



l<jj<2 

-Xnit) I dt(u'-"-^^\7U^''A-dtV'^"Ux 
+ ^3, (7.29) 

1^3 I < -^{wnitf +Wn{t) ■ Wn-l{t)) 
+ \Xn{t) - Xn-l{t)\ ■ J^'^nit). 

Adding together (j7.28p . (|7.29l) and rearranging terms, we obtain 
d /I, 



where 



2 
LI 



l<i,j<2 • 

--Ei, (7.30) 



where 



l^4| < ^(wn{t)^ +Wn{t) ■ W„-l(0) 

it) 



£2 

+ lx«(0-Xn-i(i)l-77r2^"W- (7-31) 



Now denote 



Wn{t)^l\\u'^''\t)\\l,^ + lj\{u'^-\t)y\\l. 

+ ^r<"ni)llij + ^ll(^^^"ni))'llij 

+ ^Xn(i) / I V[/(")pdx (7.32) 

+ Xn{t) E / 5,,v(""^V;")9,.C/(")dx (7.33) 

+ Xn{t) E / vf-^)9,^y/")9,,f/(")dx. (7.34) 



l<i,j<2 • 
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By using (|6.39p . (|6.2|) - ()6.4|) . it is clear that for e sufficiently small (depending on 
an absolute constant), we have 

\CM\ + \CM\ + IdLSlI 



<^\\u^''Hmh + l\\iu^''Ht)y\ 



2 

LZ 



Therefore for some absolute constant C > 0, wc have 

1 



C 



Wnit) < Wnit) < CW„{t). 



The inequalities (TTSO)) and (fTSTj) now give us 



<—{Wn{t)+Wn-,{t)) 



+ IXn(i)-Xn-i(i)| • Vt>0, 



where we have canceled a factor of Wn {t) on both sides of the inequality. 
We then inductively assume 



Therefore 



^^„.rW 



{t) 2 {t) 

£2 t t . 



< 



where in the last inequality we have used the fact that the function |x(^)— x( 2n-i )| 
is localized to the region t ^ 2". By (|6.39p . we have the a priori estimate 

We need the following Gronwall lemma 

Lemma 7.1 (Infinite time Gronwall, version 2). Let B(t) be a nonnegative smooth 
function on [l,oo) satisfying 

5 



< -B{t) , 

- t t^+<^ 



and 



B{t) < 



K 



vt > 1. 



Here < e < 1, cr > 0, 5 > 0, A' > are constants. Then for e < a. 



B{t)< 



cr - e t'' 



•-, Vi>l, 



(7.35) 
(7.36) 
(7.37) 
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where Ci{a) is a constant depending only on a. 



Proof. The proof of Lemma FTT] is similar to that of Lcmma lOl We omit the details 
here. □ 



By using Lemma 1 7. II and taking e sufficiently small, we get 

Sn+l 



Wn+lit) < 



where Sn+i ~ ■ {Sn +2 "i), where < < 1 is some constant. It is obvious that 
Sn < const-a^^ for some < a < 1. Hence the sequence ((t) (t) 2 v("+i) (t)) 

is Cauchy and converges to a nontrivial limiting function in H^. By interpo- 
lation inequalities {u^'^~^^\t),'v^'^~^^\t)) converges strongly in H"^^^ for every t. 
For the highest norm H"^, upon passing to a subsequence if necessary, we have 
(?2("+^)(i), v("+^)(i)) tends to the limiting solution (u(i),v(<)) weakly for each t. 
Hence 

\\m\H'^ + MmH^ + \\imy\\Hr^ 

+ ll(vW)'lk- 
We have obtained the desired classical solution. 
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